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Abstract 

Hexagonal circle patterns with constant intersection angles mimicking holomorphic maps z c 
and log(z) are studied. It is shown that the corresponding circle patterns are immersed 
and described by special separatrix solutions of discrete Painleve and Riccati equations. The 
general solution of the Riccati equation is expressed in terms of the hypergeometric function. 
Global properties of these solutions, as well as of the discrete z c and log(z), are established. 

1 Introduction. Hexagonal circle patterns and z c 

The theory of circle patterns is a rich fascinating area having its origin in the classical theory 
of circle packings. Its fast development in recent years is caused by the mutual influence and 
interplay of ideas and concepts from discrete geometry, complex analysis and the theory of 
integrable systems. 

The progress in this area was initiated by Thurston's idea |23|.|16| of approximating the 
Riemann mapping by circle packings. Classical circle packings consisting of disjoint open disks 
were later generalized to circle patterns where the disks may overlap (see for example |14j). 
Different underlying combinatorics were considered. Circle patterns with the combinatorics of 
the square grid were introduced in hexagonal circle patterns were studied in [7j and |S]. 

The striking analogy between circle patterns and the classical analytic function theory is 
underlined by such facts as the uniformization theorem concerning circle packing realizations of 
cell complexes with prescribed combinatorics 5 , a discrete maximum principle and Schwarz's 
lemma fF3\, rigidity properties ^H^EU and a discrete Dirichlet principle [2T] . 

The convergence of discrete conformal maps represented by circle packings was proven in 
|2U| . For prescribed regular combinatorics this result was refined. C°°-convergence for hexagonal 
packings is shown in QJ- The uniform convergence for circle patterns with the combinatorics 
of the square grid and orthogonal neighbouring circles was established in |21| . 

The approximation issue naturally leads to the question about analogs to standard holo- 
morphic functions. Computer experiments give evidence for their existence |12|.|15|. however 
not very much is known. For circle packings with hexagonal combinatorics the only explicitly 
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described examples are Doyle spirals ^^>0) which are discrete analogs of exponential maps, 
and conformally symmetric packings, which are analogs of a quotient of Airy functions For 
patterns with overlapping circles more explicit examples are known: discrete versions of exp(z), 
erf (z) |U, z c , log(z) are constructed for patterns with underlying combinatorics of the square 
grid; z c , log(z) are also described for hexagonal patterns [Zj, [Hi- 
lt turned out that an effective approach to the description of circle patterns is given by the 
theory of integrable systems (see (2] , [Jj , |Hj ) • For example, Schramm's circle patterns are governed 
by a difference equation which is the stationary Hirota equation (see [2^ ) • This approach proved 
to be especially useful for the construction of discrete z c and log(z) in P , , |H1 , with the aid 
of some isomonodromy problem. Another connection with the theory of discrete integrable 
equations was revealed in Q],[2],|3|: embedded circle patterns are described by special solutions 
of discrete Painleve II and discrete Riccati equations. 

In this paper we carry the results of for square grid combinatorics over to hexagonal circle 
patterns with constant intersection angles introduced in jSj. 

Hexagonal combinatorics are obtained on a sub-lattice of Z 3 as follows: consider the subset 

H = {{k, l,m) £ Z 3 : \k + I + m\ < 1} 

and join by edges those vertices of H whose (k, I, m)-labels differ by 1 only in one component. The 
obtained quadrilateral lattice QL has two types of vertices: for k + 1 + m = the corresponding 
vertices have 6 adjacent edges, while the vertices with k + I + m = ±1 have only 3. Suppose 
that the vertices with 6 neighbors correspond to centers of circles in the complex plane C and 
the vertices with 3 neighbors correspond to intersection points of circles with the centers in 
neighboring vertices. Thus we obtain a circle pattern with hexagonal combinatorics. 

Circle patterns where the intersection angles are constant for each of 3 types of (quadrilat- 
eral) faces (see FigQ were introduced in [S]. A special case of such circle patterns mimicking 
holomorphic map z c and log(z) is given by the restriction to an .ff-sublattice of a special isomon- 
odromic solution of some integrable system on the lattice Z 3 . Equations for the field variable 
z : Z 3 — » C of this system are: 

<l{ z k,l,mi z k,l+l,mi z k— l,Z+l,m> z k— l,l,m) =e 2lai , 

q(Zk,l,mi Zk,l,m-li Zk,l+l,m-U z k,l+l,m) = e~ 2la2 , (1) 
z k,l,m, z k+l,l,m> z k+l,l,m-li z k,l,m—l) =e > 

where on > satisfy a\ + «2 + «3 = tt and 

(Zl ~ Z2)( z 3 ~ Z A ) 



q(z 1 ,z 2 ,z 3 ,z A ) 



(Z2 ~ Z 3 )(Z4 ~ Zl) 



is the cross-ratio of elementary quadrilaterals of the image of Z 3 . Equations Q mean that the 
cross-ratios of images of faces of elementary cubes are constant for each type of faces, while the 
restriction at\ + a 2 + a 3 = tt ensures their compatibility. 

The isomonodromic problem for this system (see section [21 for the details, where we present 
the necessary results from (Hj ) specifies the non-autonomous constraint 

)( z k,l,m ~ Z k -l,l,m) , 

cz k ^ m = 2k h 

Zk+l,l,m Z k —\^l^ m 
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2| i. z k,l+l,m z k,l,m)\ z k,l,m %k,l— l,m) ^ 
2fc,Z+l,m %k,l— l,m 

n (^fcXm+l Zk.l.m) { z k,l,rn z k,Lrn— l) 
ZTO , 

z k,l,m+l z k,l,m— 1 

which is compatible with (^Q) (this constraint in the two-dimensional case with c = 1 first 
appeared in ^H)- In particular, a solution to Q,© in the subset 

Q = {(k, I, m) G 1?\ k > 0, I > 0, m < 0} (3) 

is uniquely determined by its values 

^1,0,0) z 0,l,0, z o,o,-i- 

Indeed, the constraint (J2J) gives zo,o,o = and defines z along the coordinate axis (n, 0,0), 
(0, n, 0), (0,0, — n). Then all other with (k,l,m) £ Q are calculated through the cross- 

ratios (JX|). 

Proposition 1 /<^/ T/ie solution z : Q — ► C 0/ i/ie system £|J), |]|] mi/i i/ie initial data 

2:1,0,0 = 1, 20,1,0 = e i0 , 00,0,-1 = e # (4) 

determines a circle pattern. For all (k,l,m) £ Q with even k + l + m the points %±i,z,m> z& j±i,m> 
Zklm±i li e on a circle with the center Zki t7n , i.e. all elementary quadrilaterals of the Q-image 
are of kite form. 

Moreover, equations (0) (see lemma ^ in section I3J) ensure that for the points Zk,i, m with 
k + I + m = ±1 where 3 circles meet intersection angles are Oi or 7r — a,, i = 1, 2, 3 (see Fig^ 
where the isotropic case on = ir/3 of regular and Z 3//2 -pattern are shown). 

According to proposition ^ the discrete map Zk,i, m , restricted on H, defines a circle pattern 
with circle centers Zk i m f° r k + l + m = 0, each circle intersecting 6 neighboring circles. At each 
intersection points three circles meet. 

However, for most initial data cj), i[> E R, the behavior of the obtained circle pattern is quite 
irregular: inner parts of different elementary quadrilaterals intersect and circles overlap. Define 

Q H = Qr\H. 

Definition 1 Jjjjj The hexagonal circle pattern Z c , < c < 2 with intersection angles 01,02,0:3, 
cti > 0, a\ + 02 + Q?3 = vr is the solution z : Q —* C of (0) subject to and with the initial data 

21,0,0 = 1, 20,1,0 = e ic ^ +a *\ 00,0,-1 = e ica * (5) 

restricted to Qh- 



Definition 2 A discrete map z : Qh — * C is called an immersion, if inner parts of adjacent 
elementary quadrilaterals are disjoint. 

The main result of this paper is the following theorem. 
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Figure 1: Hexagonal circle patterns as a discrete conformal map. 



Theorem 1 The hexagonal Z c with constant positive intersection angles and < c < 2 is an 
immersion. 

The proof of this property follows from an analysis of the geometrical properties of the 
corresponding circle patterns and analytical properties of the corresponding discrete Painleve 
and Riccati equations. 

The crucial step is to consider equations for the radii of the studied circle patterns in the 
whole Q-sublattice with even k + I + m. In section |31 these equations are derived and the 
geometrical property of immersedness is reformulated as the positivity of the solution to these 
equations. Using discrete Painleve and Riccati equations in section 0] we present the proof of 
the existence of a positive solution and thus complete the proof of immersedness. In section |(3 
we discuss possible generalizations and corollaries of the obtained results. In particular, circle 
patterns Z 2 and Log with both square grid and hexagonal combinatorics are considered. It is 
also proved that they are immersions. 

2 Discrete Z c via a monodromy problem 

Equations (^Q) have the Lax representation [S]: 

$k+i,l,m(.v) = ^ (1) (e,/i)^fc,«, m (^), 

**,j+l, TO (A0 = £ (2) (e,/i)$ fc ,i, m (//), (6) 
*JW,m+i(/i) = L (3) (e,/i)$ fciijm (/x), 

where \i is the spectral parameter and : Z 3 — ► GL(2, C) is the wave function. The matrices 
£(n) are d ennec l on the edges e = (p ut,Pin) of 1? connecting two neighboring vertices and 
oriented in the direction of increasing k + Z + m: 

^ (n W) = (, L Zm ~ z ° ut ), (?) 

V Zin Zout ' 

with parameters A n fixed for each type of edges. The zero-curvature condition on the faces 
of elementary cubes of Z 3 is equivalent to equations (^Q) with A n = e l5n for properly chosen 
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S n . Indeed, each elementary quadrilateral of Z 3 has two consecutive positively oriented pairs of 
edges ei,e2 and e3,e^. Then the compatibility condition 

L (m)( e2 ) L (« 2 ) (ei) = L ^Xe 4 )L^\ez) 

is exactly one of the equations Q . This Lax representation is a generalization of the one found 
in for the square lattice. 

A solution z : 1? — > C of equations Q ia called isomonodromic if there exists a wave function 
&(fj>) : Z 3 — > GL(2,C) satisfying © and the following linear differential equation in /j,: 

-j-$k,i,m(n) = A klm (p)$ klm (fi), (8) 

where Aki m (fi) are some 2x2 matrices meromorphic in /j, with the order and position of their 
poles being independent of k, I, m. 

The simplest non-trivial isomonodromic solutions satisfy the constraint: 

,2 i / ; \ ( z k+l,l,m ~ z k,l,m 

bz k,l,m + CZ k ,l,m + d=2(k- ax) h 

z k+l,l,m z k—l,l,m 

/; \ ( z k,l+l,rn z k,l,m)( z k,Lrn z k,l— l.m) . 

2(1 - a 2 ) h (9) 

z k,l+l,m z k,l— l,m 

n/ \ ( z k,l,m+l z k,l,m)(. z k,l,m z k,l'm—l) 

2{m — as) — — — — . 

z k,l,m+l z k,l,m— 1 

Theorem 2 Let z : Z 3 —> C be an isomonodromic solution to 0j with the matrix A kj i jTn in 
of the form 

r 3 B {n) 

A h>l ,M = ^ + J2^ h ^- (10) 

u;ito n-independent matrices C k ^^ m , ^ki m anc ^ normalized by tr ylo,o,o(/") = 0. Then these 
matrices have the following form: 

If -bz k j >m - c/2 bz 2 klm + cz k ^ hm + d 



^kjljin 

2 V bz k , m + c/2 



k,l,m 




- a\ 1 


z k-\-l,l,m 


z k,l,m 


{ z k+l,l,m 


z k,l,rn){ z k,l,rn 


^fc — 1,2,771 


z k+l,l,m 


z k—l,l,m * 








z k,l,m z k— 1,2,771 




R (2) 


I - 


- a 2 l 


z k,l J r l,m 


z k,l,m 


{ z k,l+l,m 


z k,l,m){ z k,l,rn 


£fc, 2 — 1,771 


z k,l+l,m 


z k,l— l,m ' 








z k,l,m z k,l— l,m 




B (3) _ 
fc,2,m 


m 


- a 3 / 


z k,l,m+l 


z k,l,m 


{ z k,l,rn+l 


z k,l,rn){ z k,l,m 


z k,l,m— 1 


z k,l,m+l 


z k,l,rn~ 1 ' 


1 






z k,l,m ^-fc, i,m—l 




2fc,2,m 


satisfies 


P. 
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Conversely, any solution z : 1? — > C to toe system is isomonodromic with Ak,i >m (n) 

given by the formulas above. 

The special case 6 = ai = 02 = 03 = with shift z ^ z — d/c implies ©. 
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3 Euclidean description of hexagonal circle patterns 

In this section we describe the circle pattern z c in terms of the radii of the circles. Such 
characterization proved to be quite useful for the circle patterns with combinatorics of the 
square grid In what follows, we say that the triangle (21, z%, Z3) has positive (negative) 

orientation if 



23 ~ zi 

Z2 ~ Z\ 



Z3 ~ Z\ 



Z2 ~ Z\ 



e i<f> with < 4> < vr (-7T < <$> < 0). 



Lemma 1 Let q(z\, 22, Z3, Z4) = e , < a < ir. 

• If \z\ — zi\ = \z\ — 24 1 and the triangle (21, Z2, 24) has positive orientation then \z3 — z%\ = 
\ z 3 — 24 1 and the angle between [#1,22] and [22,23] is (ir — a). 

• If \z\ — 22I = \z\ — 24) and the triangle (21, 22, 24) has negative orientation then (23 — 22! = 
1 23 — Z4I and the angle between [21,22] and [22,23] is a. 

• If the angle between [21,22] and [21,24] is a and the triangle (21,22,24) has positive orien- 
tation then \z3 — z<i\ = \z\ — 22I and \z3 — 24 1 = I24 — z\\. 

• If the angle between [21,22] and [21,24] is (it — a) and the triangle (21,22,24) has negative 
orientation then \ 23 — 22 1 = 1 21 — 22 1 and \ 23 — 24 1 = 1 24 — 21 1 . 

Lemma^and Proposition^imply that each elementary quadrilateral of the studied circle pattern 
has one of the forms enumerated in the lemma. 

Proposition Q allows us to introduce the radius function 

t{k,l,m) = \Zk,l,m — z k±l,l,m\ = \zk,l±l,m ~ Zk,l,m\ = \zk,l,m ~ Zk,l,m±l\, (H) 

where (k, I, m) belongs to the sublattice of Q with even k + / + m and (K, L, M) label this 
sublattice: 

k + l + m k + l + m k + l + m 
K = k , L = l , M = m . (12) 

The function r is defined on the sublattice 

Q = {(K,L,M) G Z 3 |L + M < 0, M + K < 0, K + L > 0)} 
corresponding to Q. Consider this function on 

Q H = {(K,L,M) G 7?\K > 0, L>0, M < 0, K + L + M = 0,+!}. 



Theorem 3 Let the solution z : Qh — ► C of the system J7|).p|) with initial data ^ be an 
immersion. Then function v{k,l,m) : Qh — * K+, defined by ill]), satisfies the following equations: 

(n + r2)(r 2 — r2r3 + r(r3 — r%) cos on) + (js + r2){r 2 — r2r\ + r[r\ — r%) cos a«) = (13) 

on the patterns of type I and II as in Fig\^ with i = 3 and i = 2 respectively, 

(L+M+l) + (M+K+l)— + (K+L+l) = C - 1 (14) 

r 4 + n r 6 + r 3 r 2 + r 5 
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r 2 (k,o,-k-i) r 1 



r 1 9>, 

I0L + 1.-L) \(0,L,-L-1) 





'3 0- 



(K.L+1.M) 



r '4 o r 2 / 

IL.-L) (IL.-L) J r 4 



(Kf1,L,M) 
(K,L,M) / ' 1 



Pattern II 



r \JK,L,M+1) / 

r 5 b— —4r 6 
Pattern III 




Pattern IV 



Figure 2: Equation patterns. 



on the patterns of type III and 

r(r\ sin 03 + r<i sin a\ + r% sin a^) = r\r<i sin a<i + ^r^ sin Q3 + r^r\ sin a\ (15) 

on the patterns of type IV. Conversely, r(ic,z,,M) : Qh - > ^+ satisfying equations ilcl\) ,\14\) J15\) 
is the radius function an immersed hexagonal circle pattern with constant intersection angles 
(i.e. corresponding to some immersed solution z : Qh — * C of |fl|). P))_). which is determined by 
r uniquely. 

Proof: The map 2& j )m is an immersion if and only if all triangles (Zk,i m i z k+l,l,mi z k,l,m-l)i 
( z k,l ,mi z k,l,m— l) z k,l+i,m) and (^fc,i,rri) z k+i,l,mi z k,l+i,m) of elementary quadrilaterals of the map 
z klm have the same orientation (for brevity we call it the orientation of the quadrilaterals). 

Necessity: To get equation (|14|). consider the configuration of two star-like figures with centers 
at z k,l,m with k + l + m = 1 (mod 2) and at Zk+i,i, m > connected by five edges in the /c-direction as 
shown on the left part of FigEl Let r,, i = 1, 6 be the radii of the circles with the centers at 




S Zk -'i m -' z k-i.l-!. 




Figure 3: Circles. 



the vertices neighboring Zk,i, m as i n FigEl As follows from LemmaQ] the vertices Zk,i,m-> z k+l,l,m 
and Zk—i,i,m are collinear. For immersed z c , the vertex Zk,i, m lies between Zk+i t i,m an d Zk-i,i,m- 
Similar facts are true also for the I- and m-directions. Moreover, the orientations of elementary 
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quadrilaterals with the vertex Zf. i m coincides with one of the standard lattice. Lemma ^ defines 
all angles at Zk : i : m of these quadrilaterals. Equation @ at (k,l,m) gives zj.^ m : 

c V r i + r 4 r 3 + r 6 r 2 + r 5 

where e ls = (z k+ i,i, m - Zk,l,m)/ri- Lemma [T] allows one to compute Zk+i,i,m-l, z k+i,i+i,m, 
z k+i,i,m+i and Zk+i t i-i >m using the form of quadrilaterals (they are shown in Fig. |3J). Now 
equation © at (k + 1,1, m) defines Zk+2,i,m- Condition \zk+2,i,m — z k+i,i,m\ = f\ with the labels 
JI2J) yields equation (fli]). 

For / = values z/c+i^.m, z k+2,o,m, z k+ifi,m-i and the equation for the cross-ratio with 0:3 
give the radius R with the center at Zfc+2,o,m-i- Note that for I = the term with r§ and r& 
drops out of equation (fl4"|) . Using this equation and the permutation R — > r±, r\ — > r, r 2 — > r 2 , 
r 5 r 3i one gets equation (fT3)l with i = 3. The equation for pattern II is derived similarly. 

To derive 1)15 Jl. consider the figure on the right part of FigE] where k + I + m = 1 (mod 2) 
and 7*1, r 2 , r 3 and r are the radii of the circles with the centers at Zfz+ij^, z+i m-1) z k,l+i,m 
and Zk t i t m-i, respectively. Elementary geometrical considerations and Lemma ^ applied to the 
forms of the shown quadrilaterals give equation ()15j) . 

Remark. Equation 1)15)) is derived for r = r{K,L,M), r\ = t{k,l,m-i), r 2 = r{K-i,L,u), r% = 
r(i{,L-i,M+i). However it holds true also for n = r(n,L,M+i), r 2 = v(k+i,l,m), r% = r(K,L+i,M+i) 
since it gives the radius of the circle through the three intersection points of the circles with 
radii r%, r 2 , r% intersecting at prescribed angles as shown in the right part of Fig|21 Later, we 
refer to this equation also for this pattern. 

Sufficiency: Now let t(k,l,m) : Qh — > R+ be some positive solution to (|13 )l .(|14j ) . 1)15 )1 . We can 
re-scale it so that r(0, 0,0) = 1. Starting with r(l,0, — 1) and r(0, 1, — 1) one can compute r 
everywhere in Qh'- r in a "black" vertex (see FigE} is computed from 1)14)1. (Note that only r at 
"circled" vertices is used: so to compute one needs only r(l,0, — 1) and r (0,1,-1).) The 

function r in "white" vertices on the border dQ H = {{K, 0, -K)\K £ N} U {(0, L, -L)\L G N} 
is given by (jlUjl . Finally, r in "white" vertices in Q 1 ^ 1 = Qh \ dQn is computed from ([15)1. In 
FiglU labels show the order of computing r. 

Lemma 2 Any solution t{k,l,m) : Qh R to US\) , \14\j , U5\) with < c < 2, which is posi- 
tive for inner vertices of Qh defines SOU16 satisfying f7|) in Q. Moreover^ all the tviafi- 

,mj z k+l,l,mi z k,l,m—l ), ( z k,l,mi z k,l,m—li Z k,l+X,m) O-nd ( z k,l,m, z k+l,l,m, z k,l+l,m) have posi- 
tive orientation. 

Proof of the lemma: One can place the circles with radii t(k,l,m) into the complex plane C in the 
way prescribed by the hexagonal combinatorics and the intersection angles. Taking the circle cen- 
ters and the intersection points of neighboring circles, one recovers Z\-i m for k + l + m = 0, ±1 up 
to a translation and rotation. Reversing the arguments used in the derivation of ([13 )1 .(|14 |) . 1)15 )1 . 
one observes from the forms of the quadrilaterals that equations (^Q) are satisfied. Now using 
one recovers z in the whole Q. Equation (|15|) ensures that the radii r remain positive, which im- 
plies the positive orientations of the triangles (zk,i, m , z k+i,l,m, z k,l,m-i), ( z k,i,m, z k,i,m-i, z k,l+i,m) 

and (ZkJ^m, z k+l,l,m: z k,l+l,m)- 

Consider a solution z : Q —* C of the system (|T j) .(|2 ^l with initial data (|4"|). where <j> and 
ip are chosen so that the triangles (zo,o,o)j z i,o,o, 20,0,-1) and (^0,0,0); 20,0,-1 > z o,i,o) have positive 



8 



(0,1,1 



{1,0.0) 





15 




13 



18 
16 





13 



19 
18 



15 



14 



Figure 4: Computing r in Qh- 



orientations and satisfy conditions r(l,0, — 1) = |£i,o,-i — 21,0,0 1 an d r(0, 1,— 1) = |-2o,i,— 1 — 
20,0,-1 1- The map Zk 1 m defines circle pattern due to Proposition ^ and coincides with the map 
defined by Lemma El due to the uniqueness of the solution uniqueness. Q.E.D. 

Since the cross-ratio equations and the constraint are compatible, the equations for the radii are 
also compatible. Starting with r(0, 0,0), r(l,0, — 1) and r(0, 1,— 1), one can compute r{K,L,u) 
everywhere in Q. 

Lemma 3 Let a solution t{k,l,m] : Q — > R of hl'J\) , ^T^\j , ^15\) be positive in the planes given by 
equations K + M = and L + M = then it is positive everywhere in Q. 

Proof: As follows from equation (|15|). r is positive for positive rj, i = 1, 2, 3. As r at (K, K, —K), 
(K + 1, K, —K — 1) and (K, K + 1, —K — 1) is positive, r at (K, K, K — 1) is also positive. Now 
starting from r at (K, K, -K - 1) and having r > at (N, K + 1, -K - 1) and (N, K, -K), 
one obtains positive r at (N, K, —K — 1) for < iV < K by the same reason. Similarly, r at 
(K, N, —K — 1) is positive. Thus from positive r at the planes K + M = and L + M = 0, we 
get positive r at the planes K + M = — 1 and L + M = — 1. Induction completes the proof. 

Lemma 4 Lei a solution r(n,L,M) : Q — > R of m,y\) , (I^\ ),^l,5 \) be positive in the lines parameter- 
ized by n as (n,0, — n) and (0,n, — n). Then it is positive in the border planes of Q specified by 
K + M = and L + M = 0. 

Proof: We prove this lemma for K + M = 0. For the other border plane it is proved similarly. 
Equation (JTIJ) for (K,L, -K - 1) gives 



r2 = r 5 



(2L + c)n + (2K + c)r 4 
(2if + 2 - c)n + (2L + 2 - c)r 4 



(16) 
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therefore r 2 is positive provided r l5 r 5 and r 4 are positive. For K = L it reads as 



(2K + c) 

r2 = r5 (2if + 2-c) - (17) 

It allows us to compute recursively r at (if, K, —K) starting with r at (0, 0, 0). Obviously, r > 
for (if, K, —K) if r > at (0, 0, 0). This property together with the condition r > at (n, 0, — n) 
imply the conclusion of the lemma since equation (|16|) gives r everywhere in the border plane 
of Q specified by K + M = 0. 

Lemmas |3] and ^] imply that the circle pattern z c is an immersion if r > at (N, 0, — N) and 

(0,N,-N). 



4 Proof of the main theorem. Discrete Painleve and Riccati 
equations 

In this section, we prove that all r(n, 0, — n), Vn G N are positive only for the initial data 
•2i,o,o = 1; ^oA-i = eC ° 3 - For the line r(0, n, — n) the proof is the same. Our strategy is as 
follows: first, we prove the existence of an initial value zo,o,-i such that r(n, 0, — n) > 0, Vn G N, 
finally we will show that this value is unique and is zo,o,-i = e CQ3 . 

Proposition 2 Suppose the equation 

I . i \ / 2 , \ ( x n+l + ^n/^N 2 / 2\ f x n—l + £2; n \ £ 2 — 1 . , 

(n + 1 x n - 1 — - n(l - x n /e — = (18) 

where e = e %OL ' 6 , has a unitary solution x n = e % ^ n in the sector < f5 n < 0:3. Then r(n, 0, — n), 
n > is positive. 

Proof: For -Zi 5 o.o = 1 an d unitary ^1,0.-1, the equation for the cross-ratio with 03 and @ reduce 
to (|TH)) with unitary x\ = {z n fl- n -\ — -z n) o, n )/(-2n+i,o,-n — Zn,o,n)- Note that for n = the 
term with x_i drops out of (|18j) : therefore the solution for n > is determined by a^o only. 
The condition < f3 n < 03 means that all triangles (z n fi,- n , -Zn+i,o,n> -Zn,o,-n-i) have positive 
orientation. Hence r(n, 0, —n) are all positive. Q.E.D. 

Remark. Equation ()18|) is a special discrete Painleve equation. For a more general reduction 
of cross-ratio equation see ^Hj. The case e = i, corresponding to the orthogonally intersecting 
circles, was studied in detail in pQ. Here we generalize these results to the case of arbitrary 
unitary e. Below we omit the index of a so that e = e ta . 

Theorem 4 There exists a unitary solution x n = e*^ n to i li#J) in the sector < j3 n < a. 

Proof: Equation (fTH|) allows us to represent x n+ \ as a function of n,x n _i and x n : x n+ \ = 
&(n,x n -i,x n ). &(n,u,v) maps the torus T 2 = S 1 x S 1 = {(n, v) G C : \u\ = \v\ = 1} into S 1 
and has the following properties: 

• For all n £ N it is a continuous map on Aj x Aj where Aj = {e 1 ^ : j3 S (0, a)} and Aj 
is the closure of Aj . Values of $ on the border of Aj x Aj are defined by continuity: 
<5(n,n,e) = -1, ®(n,u, 1) = -e. 



10 



• For (u,v) G Aj x Aj one has $(n,u, v) G Aj U A// U A/y, where A/j = {e i(3 : /? G (a,7r]} 
and Ajv = {e 1 ^ : (3 G [a — vr,0)}. I.e., x cannot jump in one step from Aj into Ajjj = 
{e^ :/3G (-vr,a-7r)}. 

Let xq = e ll3 ° . Then x n = x n {(5o). Define S n = {[3q : »fc(/3o) G Aj V < < n}. Then 5 n is a 
closed set since <3? is continuous on Aj x Aj. As a closed subset of a segment it is a collection of 
disjoint segments S l n . 

Lemma 5 There exists sequence {<Sn } such that: 

• Sn is mapped by x n (Po) onto Aj, 



ql(n+l) ql(n) 
°n+l C 



The lemma is proved by induction. For n = it is trivial. Suppose it holds for n. As Sn"^ 
is mapped by x n ((3o) onto Aj, continuity considerations and $(n,u, e) = —1, u, 1) = — e 
imply: x n+ i{f3o) maps S l n onto Aj U A// U A/y and at least one of the segments S l n+1 C Sn^ 
is mapped into Aj. This proves the lemma. 

Since the segments of {Sn } constructed in lemma El are nonempty, there exists (3o G S n for 
all n > 0. For this /?o, the value x n ((3o) is not on the border of Aq since then x n+ i(/3o) would 
jump out of Aj. Q.E.D.. 

Let r n and R n be the radii of the circles of the circle patterns defined by Zk,i, m with the 
centers at Z2n,o,o and Z2n+l,o,-l respectively. Constraint @ gives 

2n + c 

r n+l — —/ r ' n 
2{n + 1) — c 

which is exactly formula l|17j). From elementary geometric considerations (see Fig. EJ) one gets 




Figure 5: Circles on the border. 



r n+ i - i? n cos a 

-Kn+l — "Pi r n+l 

iin - r n+ i cos a 
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(recall that a = Q3). Define 



R n 2n + c 

Pn = , 9n{c) 



2(n + 1) - c 

and denote t = cos a. Now, the equation for the radii R, r takes the form: 

9n{c)-t Pn 

Pn+i = —FT- ( 19 ) 

Pn - tg n {c) 

Remark. Equation (|19j) can be seen as a discrete version of a Riccati equation. This is motivated 
by the following properties: 

• the cross-ratio of each four-tuple of its solutions is constant as p n +i is a Mobius transform 
of p n , 

• the general solution is expressed in terms of solutions of some linear equation: the standard 
Ansatz 

Pn = — +tg n {c) (20) 
Vn 

transforms (|19|) into 

Vn+2 + t(g n+1 (c) + l)y n +i + {t 2 - l)g n (c)y n = 0. (21) 

As follows from Theorem 01 Proposition [2] and Lemma |3J equation ()19|) has a positive solution. 
One may conjecture that there is only one initial value po such that p n > 0, Vn G N from 
the consideration of the asymptotics. Indeed, g n (c) — > 1 as n — > 00, and the general solution 
of equation (f2*T|) with limit values of coefficients is y n = c\{— l) n (l + t) n + 02(1 — t) n . Thus 
= ^j^- + tg n (c) — > —1 for ci / 0. However c\,C2 define only the asymptotics of a solution. 
To relate it to the initial value po is a more difficult problem. Fortunately, it is possible to find 
the general solution to (|2*Tj) . 

Proposition 3 The general solution to \21\) is 



T ( n +k) i >n+l- C /2 p ,3-C C-l 1 

-(ciAi F(— -n,zi)+ (22) 



w " r(n + l-f) v 1 1 v 2 ' 2 '2 

. n+i-c^n^ - C C-l 1 
+ C 2A 2 7 F(^— ,- -n,z 2 ) 

where Ai = — t — 1, A 2 = 1 — t, Zi = (t — l)/2, z 2 = — (1 + i)/2 a^d is i/ie hypergeometric 
function. 

Proof: The solution was found by a slightly modified symbolic method (see 10 for the method de- 
scription and [Sj for the detail). Here, F(a, b, c, z) denotes the standard hypergeometric function 
which is a solution of the hypergeometric equation 

z(l - z)F zz + [c - (a + b + l)z]F z - abF = (23) 
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holomorphic at z = 0. Due to linearity, the general solution of (|21|) is given by a superposition 
of any two linearly independent solutions. Direct computation with the series representation of 
the hypergeometric function 



7-1 7-1 ( , 7"1\ , , , (l-^)( 



F(l - , , 1 - x + , *) = 1 + z- 2 /v i ' + ...+ (24) 

I ° (l _ ( x + 2=1)) T k ; 



(l_2=l)( 2 -2=l)...(fc_2=l) (2=i)(i + 2=i)...(fc_i 



7— 1 - 

~2~ < 



: i _ . . . (fc - ( X + 2=1)) 



+ . 



2 

shows that each summand in 1)22(1 satisfies equation (|21[) . To finish the proof of Proposition |*3 
one has to show that the particular solutions with c\ = 0, C2 7^ and c\ 7^ 0, C2 = are linearly 
independent. This fact follows from 

Lemma 6 As n — > 00, function h2ty) has the asymtotics 

y n ^(n + l- 1 /2)^{c 1 \T 1 - 1 ' 2 + c 2 X n 2 +1 ^ /2 ). (25) 

For n — > 00 the series representation ()24j) implies F(^p-, 5 — n > z i) — !• Stirling's formula 

T(x) ~ v^e- x i x ^ (26) 

T(n+i) 

yields the asymptotics of the factor r ^ TO+1 ii- ) • This completes the proof of the lemma and of 
Proposition ©. 

Proposition 4 A solution of the discrete Riccati equation \iy\) with a 7^ n/2 is positive for all 
n > if and only if 

sin ^ , v 

Po = -72^ • ( 27 ) 

sm i — 

Proof: For positive p n , it is necessary that c\ = 0: this follows from asymptotics (|25|) substituted 
into (j2U|) . Let us define 

(l-^i)(2^i) (jfe_2=l)...(l_2=l)(2=l)(A;-l + 2=i) 

a(z) = l + z^ \ A 2 ; + + ^ ^ 2 / V 2 / V 2 —... (28) 

2 ^'V^ 2/ ' ' ' 2 

This is the hypergeometric function F(^ip-, \ — n, z) with n = 0. A straightforward com- 
putation with series shows that 

P0 = 1 + M z + ifl£^£) (29) 

2 — 7 2 — 7 

where z = =4^. Note that as a function of 2 satisfies an ordinary differential equation of first 
order since ^-pr satisfies the Riccati equation obtained by a reduction of (|23|) . A computation 



s(z) 
sin 



7 cv 



shows that — (2 _ 2 )Q satisfies the same ordinary differential eqiation. Since both expression E3) 

sin — 



and (|27|) are equal to 1 for 2 = 0, they coincide everywhere. Q.E.D. 

Proof of Theorem^ Proposition (j3J implies that the initial xq for which (|18|) gives positive r is 
unique and implies the initial values (JSJ) for z c if aj 7^ 7r/2. For the case a = tt/2, any solution 
for (fTT?|t with po > is positive. Nevertheless, as was proved in pQ, xq is in this case also unique 
and is specified by (|27jl. Thus for all n 6 N we have r(n,0, — n) > 0, r(0,n, — n) > for the 
circle pattern z c . Lemmas H] and El complete the proof. 
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5 Hexagonal circle patterns z 2 and Log 



For c = 2, formula (|17|) gives infinite r(l, 1, — 1). The way around this difficulty is renormalization 
z — > (2 — c)z/c and a limit procedure c — ► 2 — 0, which leads to the re-normalization of initial 
data (see jSj)- As follows from (|27|). this renormalization implies: 

sin a-? . , sinao , , 
r(0,0,0 =0, r 1,0,-1 = ±, r 0,1,-1 = r 1,1,-1 ) = 1. (30) 

Proposition 5 T/ie solution to Ucl\) . \14\ ), with c = 2 and initial data 13 (J\) is positive. 

Proof: This follows from Lemmas 01 and 0] since Theorem0]is true also for the case c = 2. Indeed, 
solution x n is a continuous function of c. Therefore it has a limit value as c — > 2 — and it lies 
in the sector Aj. 

Lemma |21 implies that there exists a hexagonal circle pattern with radius function r. 
Definition 3 The hexagonal circle pattern Z 2 has a radius function specified by Proposition^ 

Equations (|13 J) .([14 |) .(|15 j) nave t ne symmetry 

1 



r 



c -► 2 - c, (31) 



which is the duality transformation (see §5]). The smooth analog / — > f* for holomorphic 
functions f(w),f*(w) is: 

df{w) df*{w) _ x 
dw dw 

Note that log* (to) = w 2 /2. The hexagonal circle pattern Log is defined [H] as a circle pattern 
dual to Z 2 . Discrete z 2 and Log are the first two images in Fig. |S1 

Theorem 5 The hexagonal circle patterns Z 2 and Log are immersions. 

Proof: For z 2 this follows from Proposition |SJ Hence the values of 1/r, where r is radius function 
for z 2 , are positive except for r(0, 0,0) = oo. Lemma [2] completes the proof. 



6 Concluding remarks 

In this section we discuss corollaries of the obtained results and possible generalizations. 



6.1 Square grid circle patterns z c and Log 

Equations (^Q) extend Z\. \ m corresponding to the hexagonal z 2 and Log from Qh into the three- 
dimensional lattice Q. The r-function of this extension satisfies equation (|15j) . Consider i m for 
the hexagonal z c and Log restricted to one of the coordinate planes, e.g. I = 0. Then Proposition 
^states that Zk o, m defines some circle pattern with combinatorics of the square grid: each circle 
has four neighboring circles intersecting it at angles and it — a^. It is natural to call it square 
grid z c (see the third image in Fig. EJ). Such circle patterns are natural generalization of those 
with orthogonal neighboring and tangent half-neighboring circles introduced and studied in |21j . 

Theorem 6 Square grid z c , < c < 2 and Log are immersions. 
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Figure 6: Hexagonal z 2 , Log and square grid z 3 / 2 . 



Proof easily follows from lemma |2J 

It is interesting to note that the square grid circle pattern z c can be obtained from hexagonal 
one by limit procedure a 3 — > +0 and by a\ — > tt — a 2 . These limit circle patterns still can be 
defined by ©, © by imposing the self-similarity condition Zk i m — fi.k—m- 

6.2 Square grid circle patterns Erf 

For square grid combinatorics and a = tt/2, Schramm £2J constructed circle patterns mimicking 
holomorphic function erf (z) = (2 /it) J e~ z dz by giving the radius function explicitly. Namely, 
let n,m label the circle centers so that the pairs of circles c(n,m), c(n + l,m) and c(n,m), 
c(n,m + 1) are orthogonal and the pairs c(n,m), c(n + l,m + 1) and c(n,m + 1), c(n + l,m) 
are tangent. Then 

r(n,m) = e nm (32) 
satisfies the equation for a radius function: 

R 2 (ri + r 2 + r 3 + r 4 ) - (r 2 r 3 r 4 + rir 3 r 4 + rir 2 r 4 + rir 2 r 3 ) = 0, (33) 

where R = r(n,m), r\ = r(n + l,m), r 2 = r(n,m + 1), r 3 = r(n — l,m), r 4 = r(n,m — 1). 
For square grid circle patterns with intersection angles a for c(n,m), c(n + l,m) and it — a for 
c(n,m), c(n,m + 1) the governing equation (|33|) becomes 

R 2 (ri + r 2 + r 3 + r 4 ) - (r 2 r 3 r 4 + rir 3 r 4 + rir 2 r 4 + rir 2 r 3 ) + 2Rcosa(rir 3 - r 2 r 4i ) = 0. (34) 

It is easy to see that 1)34(1 has the same solution 1)32(1 and it therefore defines a square grid circle 
pattern, which is a discrete Erf. A hexagonal analog of Erf is not known. 

6.3 Circle patterns with quasi-regular combinatorics. 

One can deregularize the prescribed combinatorics by a projection of Z n into a plane as follows 
(see p2 ). Consider Z™ C R n . For each coordinate vector ej = (e|, ...,ef), where = 8f define 
a unit vector ^ in C = R 2 so that for any pair of indices vectors form a basis in R 2 . 
Let f2 G R n be some 2-dimensional simply connected cell complex with vertices in Z". Choose 
some xq G O. Define the map P : — > C by the following conditions: 
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• P(>o) = P 



• if x, y are vertices of f2 and y = x + e« then -P(y) = -P(x) + £j. 
It is easy to see that P is correctly defined and unique. 

We call Q a projectable cell complex if its image uj = P(£l) is embedded, i.e. intersections of 
images of different cells of f2 do not have inner parts. Using projectable cell complexes one can 
obtain combinatorics of Penrose tilings. 

It is natural to define "discrete conformal map on u n as a discrete complex immersion function 
z on vertices of uj preserving the cross-ratios of the w-cells. The argument of z can be labeled 
by the vertices x of Q. Hence for any cell of 0, constructed on e^, e,-, the function z satisfies the 
following equation for the cross-ratios: 



where ctk,j is the angle between and £j, taken positively if (£&,£?) h as positive orientation 
and taken negatively otherwise. Now suppose that z is a solution to defined on the whole 
Z™ . We can define a discrete z c : uj — > C for projectable Q as a solution to l|35j ) .l|36 )) restricted 
to f2. Initial conditions for this solution are of the form (JSJ) so that the restrictions of z to each 
two-dimensional coordinate lattice is an immersion defining a circle pattern with prescribed 
intersection angles. This definition naturally generalizes the definition of discrete hexagonal and 
square grid z c considered above. 

We finish this section with the natural conjecture formulated in [S]. 

Conjecture. The discrete z c : uj — » C is an immersion. 

First step in proving this claim is to show that equation is compatible with the constraint 



For n = 3 this fact is proven in [S]. 
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